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Abstract 



^^ In this paper we introduce a novel Mittag— LcfHcr-typc function and study its properties in relation 

to some integro-differential operators involving Hadamard fractional derivatives or Hyper-Bessel-type 

>^^ operators. We discuss then the utility of these results to solve some integro-differential equations 

^Nj involving these operators by means of operational methods. We shovi' the advantage of our approach 

through some examples. Among these, an application to a modified Lamb-Bateman integral equation is 

presented. 
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1 Introduction 



Fractional calculus is a developing field of the applied analysis concerning methods and applications of 

' ' integro-differential equations involving fractional operators. In literature there are many different definitions 

, of fractional derivatives and integrals [14]. One of these definitions was introduced by Hadamard [7] in 1892 

►^ and, although not so frequently used in applications, it is generally mentioned in classical reference books 

ryr) on fractional calculus [14, Section 18.3]. In some recent papers by Babusci et al. [2, 1], it has been shown 

l^ by means of operational techniques that the Hadamard integral corresponds to a negative fractional power 

f^ of the operator x-£^. By exploiting the integral representation of such operator, they have shown that the 

T-H Lamb-Bateman integral equation can be rearranged as a fractional integro-diflerential equation involving 

l' Hadamard fractional derivatives of order 1/2. However, a very few works on applications of Hadamard 

(^— ^ integrals and derivatives have been carried out, mainly due to the intrinsic intractability of such operators. 

rf\ On the other hand, many papers on the properties of Laguerre derivatives and their applications in partial 

,-H differential equations (see for example Dattoli et al. [3]) are present in literature. By means of operational 

IL" methods Dattoli and Ricci [4] have found a simple way to solve analitically a class of generalized evolution 

problems involving Laguerre-type operators. The purpose of this paper is to extend these results and to 

consider some applications in order to solve integro-differential equations involving Hadamard fractional 

%^ integrals and derivatives. 

Cu First, we define a general class of functions related to fractional operators involving Hadamard derivatives 

and fractional Hyper-Bessel-type operators. Some of the main properties of this class of functions and of 

the related operators will bee analyzed and discussed. Then we will show the advantage of this approach to 

solve some integro-differential equations and present a concrete example treating a modified Lamb-Bateman 

equation. 

2 Basics concepts on Hadamard fractional calculus and Laguerre derivatives 

In order to make the paper self-contained, here we recall some definitions and properties of Hadamard 
fractional integration and differentiation. We refer to Kilbas [9] for a more detailed analysis. The first step 
is to define the Hadamard fractional integral and derivative operators. 



Definition 2.1. Let 5R(q) > 0. The Hadarnard fracUonal integral of order a, applied to the function 
f € L^[a, b], < a < b < cx), x £ [a, b], is defined as 

i"/W = j^/ (log-) fit)-. (2.1) 

Before constructing the corresponding derivative operator, we define tiie following space of functions. 

Definition 2.2. Let [a, b] be a finite interval such that — oo < a <b < oo and let AC[a, b\ be the space of 
absolutely continuous functions on [a,b]. Let us denote 5 — x^ and define the space 

AC'i[a,b] = {g: [a,b] ^ C: S^^ [g(x)] e AC[a,b]} ■ (2.2) 

Clearly AC^[a, &] = AC[a,b]. Analogously to the Riemann-Liouville fractional calculus, the Hadarnard 
fractional derivative is defined in terms of the Hadarnard fractional integral in the following way. 

Definition 2.3. Let S = x-^, 5R(a) > and n = [a] + 1, where [a] is the integer part of a. The Hadarnard 
fractional derivative of order a applied to the function f G AC^[a, b], < a < b < oo, is defined as 

oyix) = r {r-y) (x). (2.3) 

A more general Hadamard-type fractional operator has been studied by Kilbas [9]. We will stick anyway to 
the classical definition. 

It has been proved (see e.g. Kilbas [9, Theorem 4.8]) that in the space L^[a,b], < a < b < cx), 1 < p < oo, 
the Hadamard fractional derivative is the left-inverse operator to the Hadamard fractional integral, i.e. 

D'^rfix) = fix). (2.4) 

Furthermore, we recall some recent results of Babusci et al. [1] on the operator {x-^) , with a > 0. The 
authors have proved by means of operational techniques that it corresponds to the Hadamard fractional 
integral as it is shown in the following. Babusci et al. [1] have considered the following integral representation 
of the power a > of a general operator a 

' ' -"" "-' (2.5) 

•- K"-) Jo 

so that 

d^ s . (2.6) 



n~" / 


1 —sa a 

ase s 


" r(a) i 


d Y" 1 


fOO 

/ Ase-'"'' 
'() 


dx) r(a) , 



To explain it clearly we recall the properties of the dilation operator e^^^x [13]. By considering x = e 
have 



„^^;fe fl^\ _ „^m t(9\ _ tte+^\ _ tl,.^ 



e 



fix) = e'^fie") = /(e"+^) = fie^x), (2.7) 



where we used the shift action of the exponential operator, that is e do y(g) = f(^0 + A). Going back to (2.6) 
^^] ■fi^) = ^J dse-=^-/(:r)s"-^ (2.8) 
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Since (2.4) holds we immediately have that (x-^) = D". This is consistent with the interpretation of the 
Hadamard fractional derivative as a fractional power of the operator 5 [14, Section 18.3]. 
These considerations, that will be clear in the following , imply for example that the Lamb-Bateman 
equation (see Section 4) can be restated as an equation involving a Hadamard fractional derivative. 
We now recall the definition of the Laguerre derivative (see e.g. Dattoli et al. [3], Dattoli and Ricci [4]) as 
the operator Dl = -^^x-^ = -^5. Furthermore, let us denote with D„l = -^^x-^ . . . g^a^g^, n G N U {0}, 
the Laguerre operator of order n and containing n + 1 derivatives. For example Dql = 4-, Dil = Dl, and 
D2L = ^x-^x-^. For the subsequent developments we also recall the following theorem by Dattoli and 
Ricci [4, Theorem 2.2] concerning the eigenfunctions of the Laguerre operator of order n. 



Theorem 2.1. The L- exponential function of order n 

oc ^ 

^"(^) = Efc^' nGNU{0}, x-GM, (2.9) 

fc=0 

ts an eigenfunction of the operator D„l, i-e. DnLen{ax) = ae„(aa;), a G K. 

The L-exponential function is also called in literature as nL-exponentials [13] and includes the two relevant 
specific cases of the classical exponential function (n = 0) and the Oth-order Bessel-Tricomi function 

(n=l). 

Definition 2.4. Let a G (0, oo), x > Q, we define the following operator of order a. 

S"/W = j-^D-fi^)- (2-10) 

Remark 2.1. Notice that when a G (0, 1), the operator (2.10) can he written as 

^"f(^) - ^ (^^) (j'-y) (^) = Dl {J'-V) (x), (2.11) 

and thus it can be seen as a generalization of Hadamard derivatives by means of Laguerre derivatives. 

In the following section we will show some advantages of this approach in order to solve explicitly a class of 
integro-differential equations involving Hadamard and Laguerre type derivatives. 

3 Equations involving the operator 3?" 

We first introduce a rather general function which will be specialized in several ways. 

Definition 3.1 (a-Mittag-Leffler function). The a-Mittag-Leffter function Ea-v,~,{x), a; G K, a > —1, 
i^ > 0, 7 G M, is defined as 



^— (-)-E r^.H^fe + 7) - ^'-'^ 

Clearly, for a = the classical Mittag-Leffler function Eo:i,^j(x) — E^^j{x) is retrieved. The following 
theorem presents the associated Laplace transform for a G (0, cx)), a; > 0. 

Theorem 3.1. We have that 
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e''"^x^''^Ec-u,',iXx'")dx^ —Ea-i-^,^ ( —] , s > 0, A G K, a; > 0, q G (0, cx)). (3.2) 
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Proof. By direct calculation we obtain 

/■oo ^ \k pec 

/ e-'''x-'-'E^.,,,^{\x'') dx = J2 ro'+Uuk + 1) e-'''x'"'+^-^dx (3.3) 
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^V A'' r{uk + j) ^ 1 fx^ 

^r"+i(i/fc + 7) 5'''=+^ s-t " ^'"'^ Vs" 

where the inversion of the sum with the integral is permitted by Theorem 30.1 of Doetsch [5]. D 

Note that when a — >■ we obtain, as expected, the Laplace transform of the two-parameters Mittag-Leffler 
function, i.e. 
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e-'"'x'<-^E^,^{\x'') Ax = ^ -, s > lAl'/". (3.4) 

s — A 



Theorem 3.2. For a = n G NU {0}, the n- Mittag-Leffler function £„;„,i(Ax''), A G M, x > 0, ;/ > 0, is 
an eigenfunction of a fractional hyper-Bessel-type operator 

d" d" d" d" 

-—x"-— -—x"-— (3 5) 

da;" da;"" ' ' ' da;'" dx" ' 
^ , ' 

n -\- 1 derivatives 

where (F jdx^ represents in this case the Caputo fractional derivative [12]. 



Proof. The statement is readily proved by simply observing that 

d" , d" d" ^ d" .^ A''a;"'° 

da::" da;" " ' dx" da;" ^ T^+ifj/fc + 1) 

^ ^ 'k=0 ^ ' 

n + 1 derivatives 

_ d" , d" d" , d" ^ A'=a;"''= r{uk + l) 

^ da;"^ da;""'da;"^ dx" -^ r"+^(iyk + 1) Fii^k + 1 - u) 

n derivatives 

-'" - da^ 2^ V^+^{uk + 1) V^{uk + I - u) 



(3.6) 
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= AS„;i,,i(Aa;'^), A G R, a; > 0. 
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Remark 3.1. Notice that i/ a = ri G N U {0} we can write the a-Mittag-Lejfler function (3.1) as a 
generalized Wright function. Indeed 



En.;v,'y{x) — 2_^ ^^^ 
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= lip' 
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with aj,bj G C, Aj,Bj G (0,cxj), 2; G C, is i/ie generalized Wright function. 

Clearly, the n-Mittag-Leffter function admits other representations, for example as a multi-index Mittag- 

Leffier function [10]. 

Remark 3.2. It is worth noticing that for {v, 7) = (1, 1) the operator (3.5) and the function (3.1) coincide 
respectively with the operator DnL and the L-exponential function (2.9). 

By means of the following definition we specialize function (3.1) for (1^,7) = (1,1) thus obtaining a 
generalization of the L-exponential function (2.9) for fractional values of parameter n. 

Definition 3.2. The aL- exponential function ta{x), a G (— l,cx)), x G K, is defined as 



ec{x) = J2 



fc!"+i' 



(3.9) 



For the Cauchy-Hadamard theorem, the aL-exponential function is clearly convergent in K for a G (—1, 00) 
Its Laplace transform for x > 0, a G (0, cx)) reads 



/ e^°^ec,(a;) da; = - yj 
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(3.10) 



where the inversion of the sum with the integral is permitted by Theorem 30.1 of Doetsch [5]. Note that 
formula (3.10) for a — > 0, s > 1, correctly gives the Laplace transform of the classical exponential function. 
Note that, as expected, the above result is consistent with (3.2). 
Figures 1 and 2 show the aL-exponential function for some values of parameter a. 

Theorem 3.3. The aL- exponential function ec,(Aa;), a G [0,oo), a: > 0, A G K, is an eigenfunction of the 
operators 



— a;— — a;— D"+^"'' 
da; da; da: da; 



Q G [0,00), r — 1 



,•••,'", 



(3.11) 



r derivatives 



where n = [a + 1] is the integer part 0/ a + 1. Moreover, for a G (—1, 0) the aL- exponential function ea,(Aa;), 
a: > 0, A G K, is on eigenfunction of the operator x^^D"'^^ . 




Figure 1: The aL-exponential function ea{x), a; > in loglinear plot (left) and linear plot (right), for a = 
(black — classical exponential), a = 0.2 (green), a = 0.4 (red), a = 0.6 (blue), a = 0.8 (magenta), q = 1 
(cyan — Oth-order Bessel-Tricomi function). 



Proof. For the first part we can write that 



dx dx dx dx 

r derivatives 



taO^x) = -;-K-— . . . -r-X—- 

ax Ax Ax Ax 



A A x^,fefc"+i-''a:'= 



E^^ 



k\^ 



(3.12) 
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Analogously, in order to prove the second part it is sufficient to note that 

k^ 
~fc! 



^D"+\c.{\x) =x-'^22^'' t.|a+i =^e«(^a:) 



(3.13) 
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Remark 3.3. An interesting specific case of Theorem 3.3 is when r = 1. In this case the aL-exponential 
function is an eigenfunction of the operator ■§^D°' = S" (see Definition 2.4)- This can be explained also by 
recalling that 



^ r^f^ + l — T" X^~^ r^ck + 1 — r ^ rr — 1 en — r + 1 jn — a ^ r^a T^^ 

Ax Ax Ax Ax Ax Ax Ax 

r derivatives 
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In order to describe some applications we first recall the following result [4, page 490, Theorem 5.1] 
specialized for n = 1. 

Theorem 3.4. Consider the problem 



Dl,,: Six,t) = -§-^S{x,t), x>0,t>0, 
S{0,t) = s{t), 

with the analytic boundary condition s{t). The operational solution is given by 

d' 



S{x,t) = ei 



dt 
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(3.15) 



(3.16) 



The operational solution (3.16) becomes an effective solution when the series converges, and this depends 

upon the actual form of the boundary condition S{Q,t). Notice that, from a physical point of view, this 

Laguerre-heat equation is in fact a diffusion equation with a space-dependent diffusion coefficient. 

In a similar way we can solve more general fractional partial differential equations. We notice that the 

utility of operational methods to solve fractional differential equations was pointed out by different authors, 

see for example Hilfer et al. [8], Luchko and Gorenflo [11], Garra and Polito [6], Tomovski et al. [15]. See 

also the book by Yakubovich and Luchko [16] and the references therein. 

Here we present the following result for an initial value problem involving the operator S". 




Figure 2: The aL-exponential function ea{x), a; > in loglinear plot for a = (black — classical exponential), 
Q = -0.1 (red), a = -0.2 (green), a = -0.3 (blue). 



Theorem 3.5. Consider the following initial value problem. 

's?/(x,t) = e,/(x,i), f>o, 

f{x,0)^gix), 



(3.17) 



for a G (0, 1), with an analytic initial condition g{x) and where Q^ is a generic linear integro- differential 
operator with constant coefficients, acting on x, and which satisfies the semigroup property, that is Oa;0y = 
Qx+y By applying Theorem 3.3 for r = 1, a G (0, 1), the operational solution to (3.17) is given by 



"^ j-T r^r 

f{x,t) = i^{tQx)g(x) ^^:^:^g{x)- 

Again, when the operational solution (3.18) converges, it becomes the effective solution. 
Corollary 3.1. Consider the following boundary value problem. 

'Tl"w{x,t) = Etw{x,t), a e (0, 1), t > 0, a; > 0, 
w{0,t) = h{t), 



(3.18) 



(3.19) 



with an analytic boundary condition h(t) and where Ht is a generic linear integro- differential operator with 
constant coefficients, acting on t, and which satisfies the semigroup property. Analogously to Theorem 3.5, 
the operational solution to (3.19) is given by 



,{x,t) = eo,{xEt)h{t) = Yl n^T^i)- 



(3.20) 



It is important to remark that the operational solutions (3.18) and (3.20) are indeed effective solutions 
whenever the corresponding series are convergent. 

Example 3.1. Consider the following initial value problem. 

'd? u{x,t)^ ^u{x,t), aG(0,l), i>0, 
u(x, 0) = sim. 

By Theorem 3.5 the solution is given by 

,k 



(3.21) 



u{x,t) = sin a; ea(i) = sin a; \^ —r^ 



fc!"+i' 



(3.22) 



which is clearly convergent for each {x,t) G R x R"*". 



4 An application to a modified Lamb— Bateman equation 

In [1] the authors have shown that a modified Lamb-Bateman integral equation can be written as a 
fractional equation involving Hadamard derivatives of order 1/2. They have studied the equation 



j°^ dyu{e-^\)^f{x), (4.1) 



/o 
and therefore 



where u(x) is the unknown function that has to be determined and f[x) is a known continuous function. 
Recalling that, for the dilation operator E{\) — e'^^as^ we have 

E{\) u{x) = e^^'^uix) = u{e^x), (4.2) 

then equation (4.1) can consequently be written as 

' Aye-^"^Au{x) = f{x), (4.3) 

that is, an integral equation involving Hadamard derivatives of order 1/2. 

Note that, in the general case, the function / depends on u. The simplest choice is f{u{x)) = u(x). In this 

specific case we have 

u{x) = -^D^'^u{x). (4.5) 

An application of the x-derivative to both sides of the above equation leads to a formulation of the 
Lamb-Bateman equation involving the operator D". 

'D^'^u{x) = ^-^u{x). (4.6) 

By using a simple ansatz, that is u{x) = x^ , where /? is to be specified, we find by substitution 



SO that u{x) = a;*^^ is a simple exact solution if /3 = tt/4. 
In general we can consider the equation 



/o 
which can be rewritten, as before, as 



dyu(^e-^'^''x^=f{x), MG(0,c»), (4.8) 



/ dye-" ""^uix)^ f{x). (4.9) 
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By evaluating the integral we then have that 



"(-) - rVTT) (4)'^(^) = !VTI)^''^(^)- ^'-''^ 



Again, the simplest choice of f{u{x)) = u{x) leads to 



u{x) = ^^^^^^ D'^uix), M G (0, (»). (4.11) 



The related equation involving the operator S™ is 



'D''u{x)=T[^l^l)-^u{x), (4.12) 



and by choosing again u{x) = x^ we obtain that j3 = [r(^ + 1)] 



1/m 
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